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چکیده 


پیش‌بینی دبی رسوبی با دقت بالاتر, از مهم‌ترین مولفه‌های فرآیندهای هیدرولوژیکی در مدیریت منابع آب می‌باشد. به دلیل پیچیدگی پدیده انتقال 
رسوب و وجود چندین پارامتر موثر در تخمین آن» تعیین روابط حاکم مشکل می‌باشد. اين تحقیق به منظور ارائه مدل بهینه برآورد بار رسوب معلق در دو 
ایستگاه هیدرومتری متوالی یک رودخانه طبیعی انجام گرفته است. در این راستاه جهت کاهش مقدار خطا در پیش‌بینی دبی رسوبی از روش‌های 
پیش‌پردازش سری زمانی به همراه روش‌های هوشمند مبتنی بر کرنل ماشین بردار پشتیبان (5۷۷) و رگرسیون فرآیند گاوسی (0۳15)) استفاده شده 
است. دو سناریو بررسی دبی رسوب معلق برای حالت تک ایستگاهی و ارتباط بین ایستگاهی در نظر رفته شد و مدل‌های متفاوتی بر اساس مشخصات 
هیدرولیکی و ذرات رسوبی تعریف گردید و مورد ارزیابی قرار گرفت. در استفاده از روش‌های پیش‌پردازش, ابتدا روش تبدیل موجک (۷1) به‌کار رفت» 
سپس زیر سری‌های با فرکانس بالای بدست آمده از روش ۷/۲ با روش تجزیه مد تجربی («۳(۷۲1) دوباره تجزیه گردیدند. در نهایت زیر سری‌های 
تاثیرگذار به‌عنوان ورودی مدل‌های مبتنی بر کرنل استفاده شدند. نتایج حاصل از تحلیل مدل‌های تعریف‌شده, دقت بالای روش‌های تلفیقی به کار رفته 
در تحقیق را در تخمین رسوب معلق به خوبی نشان داد. نتایج نشان داد که در حالت اول مقدار معیار خطا برای مدل برتر به ترتیب از ۰,۰۳۵ و ۰,۰۳۷ 
برای روش‌های 0۳1 و ٩۷۷]‏ به ۰,۲۸ و ۰,۲۹ برای مدل‌های تلفیقی کاهش یافت. به طور کلی» مدل‌های ترکیبی دقت مدل‌سازی را بین ۲۰ تا ۲۵ 
درصد افزایش دادند. جهت ارزیابی قابلیت اطمینان مدل برتر از تحلیل عدم قطعیت مونت کارلو استفاده شد و نتایج نشان داد که مدل 0۳3 دارای 


درجه عدم اطمینان مطلوبی در مدل‌سازی است. 


واژه‌های کلیدی: پیش پردازش, تجزیه مد تجربی» رگرسیون فرآیند گاوسی» رسوب معلق 


مقدمه 

با توجه به محدودیت منابع آبی در طبیعت توزیع زمانی و مکانی 
غیریکنواخت منابع آبیء افزايش آلودگی‌ها؛ تخریب منابع طبیعی و 
همچنین افزایش جمعیت» رشد و توسعه جوامع شهری و فعالیت‌های 
کشاورزی و صنعتی» ضرورت برنامه‌ریزی و مدیریت صحیح منابع آبی 
اجتناب‌ناپذیر است. در اين میان آگاهی از میزان مواد رسوبی که توسط 
جریان حمل یا ترسیب می‌گردد یکی از مسائل مهم و کاربردی در 
مدیریت منابع آب و از اطلاعات ضروری برای اجرای پروژه‌های آبی 
است (۸). تخمین صحیح رسوبات حمل شده به‌وسیله رودخانه در 


- گروه عمران» دانش‌کده فنی و مهندسی» واحد اهر دانشسگاه آزاد اسلامی. آهر 
ایران 
زان 


طیف وسیعی از مسائل از قبیل طراحی مخازن سدهاء انتقال رسوب و 
آودگی در دهاش ناجه‌ها و مضب‌هار طارص کاناز های نایدا 
ارزیابی زیست‌محیطی پروژه‌ها و مدیریت آبخیز به‌کار گرفته شده است 
و می‌تواند موجب کنترل بهتر عواقب زیان‌بار این پدیده و اسستفاده 
بهینه از تأسیسات آبی (مانند سده کانال -.) و جلوگیری از صرف 
هزینه‌های اضافی شود. مواد رسوبی در رودخانه‌ها به دو صورت کلی 
بار معلق و بار بستر می‌باشد که محاسبه مقدار هر یک به روش خاصی 
انجام می‌پذیرد. از جمله این روش‌ها مدل‌های تجربی و اسستفاده از 
تکنیک رگرسیون است که در این روش یک يا چند منحنی بر 
داده‌های رسوب برازش داده می‌شود (۱۱). ساده‌ترین شکل معادلات 
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رگرسیونی که مستلزم تغییر شکل لگاریتمی داده‌های دبی رسوبی و 
دبی جریان است به شکل "09-200 بوده که در آن 05 دبی رسوب و 
0 دبی جریان آب و 2 و « پارامترهای رگرسیونی می‌باسند (۱۳). با 
وجود اینکه مدل‌های مختلفی برای پیش‌بینی رسوب معلق پیشنهاد 
شده است اما با توجه به ماهیت غیرخطی و ویژگی‌های متغیر زمانی و 
مکانی در سیستم گردش آبی. هیچ‌یک از مدل‌های آماری و مفهومی 
پیشنهاد شده به‌منظور الگوسازی دقیق رسوب نتوانسته‌اند به‌عنوان یک 
مدل برتر و توانا شناخته شوند و عدم قطعیت در پیش‌بینی‌ها از اعتبار 
برآوردها از طریق مدل‌های مختلف می‌کاهد. از طرفی امروزه محققین» 
با پپشرفت روش‌های هوشمند. در جستجوی راه‌هایی برای پیش‌بینی 
پارامترهای مهم هیدرولیکی و هیدرولوژیکی هستند. در دهه‌های اخیر 
روش‌های نرم برای شبیه‌سازی رفتار سیستم‌ها در مهندسی آب مورد 
استفاده قرار گرفته‌اند. این روش‌ها اکثراً الهام گرفته از پدیده‌های 
موجود در طبیعت هس‌تند و به‌عنوان ابزاری توانمند در حل مس‌ائل 
پیچیده به شمار می‌آیند (۳). به‌عنوان نمونه سیواپراگسم و لیونگ (۱۵) 
از مدل ماشین بردار پستیبان برای پیش‌بینی بارش" رواناب حوضه 
60 واقع در کشور هلند استفاده کردند. ذو و همکاران 
(۱۸) از مدل رگرسیون فرآیند گاوسی برای پیش‌بینی دبی رودخانه 
استفاده کردند. معروف‌پور و همکاران (۱۰) از روش‌های هوش 
مصنوعی برای شبیه‌سازی تبخیر-تعرق استفاده کردند. روشنگر و 
قاسم‌پور (۱۴) با استفاده از روش‌های هوشمند رگرسیون بردار پشتیبان 
و سیستم استنتاجع عصبی- فازی انطباقی مقاومت جریان را در 
رودخانه‌های مارییچی بررسی کردند. حیدری و همکاران (۷) از روش 
هوشمند مبتنی بر کرنل برای مدلسازی مشخصات پرش هیدرولیکی 
استفاده کردند. عاشقی و همکاران (۴) از روش‌های هوش مصنوعی 
برای شبیه‌سازی بار بستر استفاده کردند. 

بررسی منابع صورت گرفته. نشان از کاربرد روش هوشمند شبکه 
عصبی مصنوعی در پیش‌بینی سری‌های زمانی از جمله رسوب 
می‌باشد. این در حالی است که مطالعات چندانی در زمینه کاربرد 
روش‌های مبتنی بر کرنل مانند رگرسیون فرآیند گاوسی در تخمین 
رسوب معلق انجام نشده است. از طرفی» در سال‌های اخیر استفاده از 
روش‌های تجزیه سیگنال به‌عنوان روشی نوین و موّثر در زمینه تحلیل 
سیگنال‌ها و سری‌زمانی مورد توجه قرار گرفته است (۱). در این روش 
سیگنال اصلی به چندین زیرسری تجزیه شده و داده‌های حاصل از 
تجزیه با معادلات مناسب برازش می‌یابد و مدل به‌دست‌آمده برای 
پیش‌بینی استفاده می‌شود. از جمله این روش‌ها تبدیل موجک و روش 
جدید تجزیه مد تجربی یکپارچه کامل است (۵) همچنین» ترکیب 
پیش‌بینی‌های چندین مدل رگرسیونی با استفاده از یک الگوریتم 
یادگیری معین می‌تواند باعث بهبود نتایج خروجی روش‌های 
رگرسیونی شود. 

بدین ترتیب با توجه به پیچیده بودن پدیده رسوب و تأثیر 


پارامترهای مختلف در پیش‌بینی آن؛ بهره‌گیری از روش‌های جدید که 
بتوانند قطعیت مدل را افزایش دهند ضروری به نظر می‌رسد. لذاء 
هدف از این تحقیق بررسی کارایی روش‌های تلفیقی پیش‌پردازش 
به‌عنوان روش جدید جهت بهبود دقت پیش‌بینی دبی رسوبی است که 
می‌تواند در برنامه‌ریزی و مدیریت بهینه منابع آبی مورد اسستفاده قرار 
گیرد. با استفاده از روش‌های تجزیه سیگنال حجم محدودی از نویز 
سفید وارد سیگنال اصلی شده و با استفاده از جنبه‌های مثبت آماری 
نویز سفید که توزیع متوازن در حوزه فرکانس دارده اثر نویز متناوب از 
روند تجزیه حذف خواهد شد. در این تحقیق» از روش‌های تبدیل 
موجک گسسته (۷۷) و تجزیه مد تجربی (1۷1۳) برای پیش‌پردازش 
مدل‌ها استفاده شده است. لذء جهت بررسی کارایی روش‌های به کار 
رفته در تحقیق, مدل‌هایی بر اساس پارامترهای هیدرولیکی و دبی 
رسوبی برای دو ایستگاه هیدرومتری واقع در رودخانه هوستونیک 
تعریف شده و دو حالت تخمین بار رسوبی بر اساس داده‌های خود 
ایستگاه و ایستگاه قبلی مورد بررسی قرار گرفته است. همچنین میزان 
قابل اطمینان بودن نتایج مدل‌های به کار رفته با استفاده از روش عدم 


قطعیت مونت‌کارلو مورد بررسی قرار گرفت. 


مواد و روش‌ها 
منطقه موردمطالعه 


داده‌های مورد استفاده در این تحقیق دبی جریان و بار معلق 
روزانه‌ی دو ایستگاه هیدرومتری متوالی از رودخانه‌ی هوستونیک 
(منصماددم) واقع در ایالات متحده آمریکا در بازه زمانی تقریبا دو 
ساله می‌باشد. اولین ایستگاه نزدیک برینگتون بزرگ در ایالت 
ماساچوست و ایستگاه دوم در روستای فالز در ایالت کانکتیکات واقع 
شده و فاصله‌ی بین آن‌ها حدود ۵۰ کیلومتر است. مساحت حوضه به 
ترتیب برای ایستگاه‌ها ۲۸۳ مایل مربع و ۶۳۴ مایل مربع می‌باشد و 
مسیر جریان از ایستگاه اول به سمت ایستگاه دوم است. داده‌های 
مربوط به بار رسوب روزانه برای منطقه مطالعاتی در سایت 1505] در 
دسترس است که اطلاعات کافی در مورد دبی رسوب و دبی جربان ر 
به دست می‌دهد. داده‌های مشاهده شده‌ی دو ایستگاه در حدود تقریبا 
دو سال از ۱ آوریل ۲۰۰۰ تا ۲۰ سیتامبر ۲۰۰۱ می‌باشد. این داده‌ها 
به‌صورت روزانه و به طور پیوسته برداشسته شده و حدود ۵۴۹ داده 
می‌باشند. موقعیت ایستگاه‌ها در شکل نشان شده است. 


روش‌های هوشمند 

ضعف عملکرد روش‌های استخراج ویژگی خطی برای داده‌های 
پیجیده که به‌صورت خطی تفکیک‌پذیر نمی‌باش ند سبب پیدایش 
روش‌های غیرخطی از جمله روش‌های مبتنی بر کرنل برای این منظور 


شده است. 
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شکل ۱- موقعیت ایستگاه‌های رودخانه هوستونیک 
کصمناهای ۱۱۲۵ »نصرماعتم۲۱ عط) ۵۴ حمنای1۱۵ 11۳6 -1 ۲1۵۱۲۵ 


در روش‌های مبتنی بر کرنل ابتدا داده‌ها به فضایی با ابعاد بیشتر 
نگاشت یافته که در این فضا به‌طور خطی تفکیک‌پذیر می‌شوند. این 
روش‌ها بر این ایده استوارند که ضرب داخلی داده‌های نگاشت يافته با 
یک تابع غیرخطی را می‌توان با تابع کرنل تخمین زد که به آن حقه 
کرنل می‌گویند. از جمله روش‌های مبتنی بر کرنل می‌توان به ماشین 
بردار پشتیبان (۷۱۷) و رگرسیون فرایند گاوسی (0۳۴1) اشاره کرد. 
الگوریتم 5۷۷ جزء الگوریتم‌های تشخیص الگو دسته‌بندی می‌شود 
رای اولین ,بار #وستط رفک هشال ۱۳۵ مش فی فت ۶ 
به‌صورت آموزش (1220) و آزمون (1600) می‌باشد. مزیتی که 5۷ 
نسبت به سایر الگوريتم‌ها دارد این است که این الگوریتم مبنی بر 


احتمالات نمی‌باشد و همواره جواب یکه‌ای به همراه دارد و زمان 
کمتری را برای دستیابی به جواب بهینه به خود اختصاص می‌دهد. از 
الگوریتم ٩۷۷‏ در هرجایی که نیاز به تشخیص الگو يا دسته‌بندی 
اشیا در کلاس‌های خاص باشد می‌توان استفاده کرد. درواقع مبنای 
کاری ٩۷1۷۲‏ دسته‌بندی خطی داده‌ها است و در تقسیم خطی داده‌ها 
سعی می‌شود خطی انتخاب شود که حاشیه اطمینان بیشتری داشته 
باشد. هدف اصلی الگوریتم ٩۷۷۲‏ این است که خطی را بیابد که از 
داده‌های موجود در دو کلاس دارای بیشترین فاصه باشد و یا 
به‌عبارت‌دیگر دارای کمترین ریسک عملیاتی باشد. فرایند گاوسی نیز 
روشی مفید به‌منظور تعریف توزیع‌های مقدم برای مدل‌های 
انعطاف‌پذیر رگرسیون و طبقه‌بندی است که در آن رگرسیون یا توابع 
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اختمال.طیقه به اشتکال سانه پارامتری محنود نف ده پاش سید قرایند 
گاوسی محموعه‌ای از متغیرهای تصادفی است که تعداد دلخواه 
محدودی از آن‌ها با توزیع‌های گاوسی ادغام شده‌اند. فرایند گاوسی 
به‌طور کامل توسط تابع میانگین (00: و تابع کوواریانس (۶ ,1:6 
آن تعیین می‌شود. اين فراینده تعمیمی طبیعی از توزیع گاوسی است 
که میانگین و کوواریانس آن به ترتیب یک بردار و یک ماتریس است. 
مدل‌های رگرسیون فرایند گاوسی بر مینای این فرض هستند که 
شتا هدات کی بای شام اطعا دوز هلگ با شا ی از 
جذابیت‌های فرایند گاوسی تنوع بسیار توابع کوواریانس آن است که به 
انخاد تولیعی با درجات:متلت.همواری با انوا مطلفه بت اغتازهای 
پیوستی منجر می‌شود و به محقق این امکان را می‌دهد که از میان 
آن‌ها گزینش مناسبی انجام دهد. اين مدل‌ها می‌توانند توزیع‌هایی را 
بین توابعی با تعداد یک یا چندین متفیر ورودی مشخص کنند. 
هنگامی که چنین تابعی میانگین پاسخ را در یک مدل رگرسیونی با 
خطاهای گاوسی تعریف می‌کند. می‌توان از محاسبات ماتریس برای 
استنتاج بهره‌گیری کرد؛ این مسئله برای محموعه داده‌هایی با بیش از 
هزار نمونه امکان‌پذیر است. طراحی روش‌های مبتنی بر کرنل شامل 
استفاده از مفهوم تابع کرنل است. در واقع با یک تبدیل غیرخطی از 
فضای ورودی به فضای خصیصه با ابعاد بیشتر (حتی نامتناهی) 
می‌توان مسائل را به‌صورت خطی تفکیک‌پذیر ساخت. از مهم‌ترین 
توابع کرنل می‌توان به کرنل خطی» چندجمله‌ای ساده. چندجمله‌ای 
نرمال‌ ده تابع شعاع محور, تابع کرنل پیرسون و سیگموئید اشاره 
کرد. 


تبدیل موجک 

موجک (/۷۷2۷۵16) دسته‌ای از توابع ریاضی است که برای تجزیه 
سیگنال پیوسته به ملفه‌های فرکانسی آن بکار می‌رود که سطح 
تفکیک هر موّلفه برابر با مقیاس آن است. تبدیل موجک تجزیه یک 
تبع بر مبنای توایع موجک است. موجک‌ها (که به‌عنوان موجک‌های 
مادر شناخته می‌شوند) نمونه‌های انتقال‌یافته و مقیاس شده یک تابع 
[موجک مادر) با طول متناهی و نوسانی شدیدا میرا هستند. تابع موجک» 
تابعی است که دو ویژگی مهم نوسانی بودن و کوتاه‌مدت بودن را دارد. 
()/۱ تابعی موجک است اگر و فقط اگر تبدیل فوریه آن (۱۳)۵ شرط 
زير را ارضا کند (0). 
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)۱( مب > 6 سس 
این شرط با عنوان شرط پذیرفتگی برای موجک (۱6 شناخته 


می‌شود ()۷۷ تابع موجک مادر است (رابطه ۲) که ضرایب ‏ و ۵ در 
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توابع موجک دارای انواع بسیاری هستند که مهم‌ترین و 
پرکاربردترین آن‌ها شامل تابع موجک 11227 2۷676۲ نقط صدءن۷ 
٩۱6۶ 6‏ فمنطمعطانیدنآ اعزم) اسست. تبدیل موجک 
پیوسته به‌عنوان روش جایگزین تبدیل فوریه زمان-کوتاه ارائه گردید و 
هدف آن, فائق آمدن بر مشکلات مربوط به تفکیک‌پذیری در تبدیل 
فوریه زمان -کوتاه است. در عمل اکثر سری‌های زمانی هیدرولوژیکی 
به‌صورت گسسته در زمان‌هایی مشخص اندازه‌گیری می‌شوند. بنابراین 
کاربرد موجک گسسته می‌تواند تناسب بهتری با این نوع سری زمانی 
داشته باشد. تبدیل موجک گسسته در عمل همان تبدیل موجک پیوسته 
است (نسخه گسسته شده آن) که در آن یک سری موجک از تبدیل 
موجک پیوسته نمونه گرفته است و سیستم‌های کامپیوتری قابلیت 
محاسبه آن را با هزینه زمانی کمتر دارا است. اصول تبدیل موجک 
گسسته به روشی تحت عنوان کدینگ زیرباند برمی‌گردد که در سال 
۶ سنگ بنای اولیه آن گذارده شد و در آن نوعی توصیف زمان- 
مقیاس از سیگنال گسسته با استفاده از فیلترهای دیجیتال ارائه می‌گردد. 
تبدیل. موجک» حاصل شیاهت‌سنجی, (کورولیشن! بین. محتوای 
فرکانسی (مقیاسی) سیگنال و تابع موجک در مقیاس‌های مختلف 
است. تجزیه کردن سری‌های زمانی با استفاده از فیلترهای مشخص 
(موجک و توابع مقیاس گذاری شده) دو نوع زیرسری را به وجود 
می‌آورند؛ بردارهای تقریبی یا باقی‌مانده و بردارهای جزییات )٩(‏ این 
ضرایب حاصل از در هم پیچیدن سری زمانی اصلی با فیلتر پایین -گذر 


ای وفناید بااگوز مات تیم مک شتسه نع تیال موعک 
قابلیت تجزیه سری زمانی به چندین زیرسری زمانی با مقیاس‌های 
مختلف را دارد و با بررسی زیرسری‌های زمانی منتجه از سری زمانی 
کلی» رفتار کوچک مقیاس و بزرگ مقیاس یک فرآیند هیدرولوژیکی 
را آنالیز می‌کند. روابطی که بیشتر برای تخمین سطح تجزیه مورد 
استفاده قرار می‌گیرد به‌صورت زیر است که در آن .1 سطح تجزیه. ۷ 
تعداد نقاط داده در سری زمانی و 0 برابر با زمان محو شونده برای 
موجک مادر انتخاب شده است. 
از 
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روش تجزیه مد تجربی ۲ ۲۷۲۷ 
روش تجزیه مد تجربی روشی برای تجزیه سیگنال‌های گوناگون 
است که این عمل در فرایندی تحت عنوان غربال کردن صورت 


سیدمهدی اقبیان. تخمین بار معلق رسوبی با استفاده از روش‌های هوشمند تلفیقی با در نظر گرفتن عدم قطعیت مدل ۴۷۹ 


می‌گیرد. طی اين فرایند سیگنال اصلی به تعدادی مولفه با محتوای 
بسامدی متفاوت تجزیه می‌شود. روش 15۱۷/9 طبق رابطه (۴) سیگنال 
اصلی () را به تعدادی مد ذاتی ("1]۳) تجزیه می‌کند (۶). 
1 
)۳ 00 10 2۰:00 < (۲ 
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که (:۲ همان مولفه باقی‌مانده. بعد از 0 تعداد 11۳5 و (ه 
است. یک داده ممکن است در یک زمان دارای چندین مد ذاتی باشد. 
این مدهای نوسانی تابع‌های مد ذاتی (171۳) نام دارند و دارای دو 
شرایط هستند: ۱) در کل داده‌هاء تعداد نقاط فرین (اکسترمم) و نقاط 
نقطه میانگین پوش برازش داده شده بر نقاط بیشینه محلی و پوش 
برازش داده شده بر نقاط کمینه محلی باید صفر باشد. به علّت وجود 
تناوب و نویز در سیگنال‌ها, در برخی موارد به علّت اختلاط مدها؛ 
توزیع حوزه زمان- فرکانس دچار انقطاع می‌شسود و عملکرد ۲11 
دچار نقض می‌گردد (چون میانگین ۷۳آها حالتی نامشخص به خود 
می‌گیرد). جهت رفع اين ایراده وو و هوانگ (۲) روشی متفاوت با عنوان 
روش تجزیه مد تجربی یکپارچه (212۷117) را پیشنهاد کردند. در روند 
تجزیه‌سازی آن» حجم محدودی از نویز سفید وارد سیگنال اصلی 
می‌شود. با استفاده از جنبه‌های مثبت آماری نویز سفید که توزیع 
متوازن در حوزه فرکانس دارد اثر نویز متناوب از روند تجزیه حذف 
می‌شود. 


معیارهای ارزیابی مدل‌ها 

به‌منظور ارزیابی کارایی روش‌های استفاده‌شده در تحقیق, از سه 
پارامتر آماری استفاده گردید که عبارت‌اند از ضریب همبستگی بین 
مقادیر مشاهداتی و محاسباتی (415 ضریب تبیین () و ريشه 
میانگین مویسات خطاها 065۳ هر چه سقذار ۱0 پرای ,نک متل 
به یک نزدیک‌تر و مقدار 1۷1517 کوچک‌تر باشد به معنی مطلوب 


بودن آن مدل است. روابط این پارامترهای آماری به‌صورت زیر 
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در روابط بالا ما: مقدار اندازه‌گیری‌شده مبا: متوسط مقدار 
اندازه‌گیری‌شده» سبا: مقدار پیش‌بینی‌شده وبا: متوسط مقدار 
پیش‌بینی‌شده و آ تعداد داده‌هاست. واردکردن داده‌ها به‌صورت خام 
باعث کاهش سرعت و دقت شبکه می‌شود لذا نرمالیزه کردن داده‌ها کمک 
شایانی به آموزش بهتر و سریع‌تر مدل می‌کند. در تحقیق کنونی داده‌ها 
مطابق رابطه (۸) نرمالیزه شدند. که در آن مسا و سم به ترتیب حداقل و 
حداکثر مقدار داده‌های مشساهده‌شده و مد داده نرمال شده می‌با د. 
محدوده‌ی نرمالیزه کردن داده‌ها در این تحقیق بین 9۰/۱ ۱ انتخاب گردید. 


)۸( | لت ]ده :ه- ورد 


صتصصر؟ - یمور 


عدم قطعیت مونت کارلو 

به صورت کلی» شبیه سازی مونت کارلو به هر تکنیکی اتلاق 
می‌شسود که از طریق نمونه‌سازی آماری» پاسخ‌های تقریبی برای 
مسائل کمی فراهم می‌کند. شبیه‌سازی مونت کارلو بیشتر برای 
توصیف روشی جهت انتشار عدم قطعیت‌های موجود در ورودی مدل 
به عدم قطعیت‌ها در خروجی مدل, به کار می‌رود. بنابراین مونت‌کارلو, 
نوعی شبیه‌سازی است که به طور صریح و کمی, عدم قطعیت را 
نمایش می‌دهد. شبیه‌سازی مونت کارلو متکی به فرآیند نمایش صریح 
عدم قطعیت با تعیین ورودی‌ها به عنوان توزیع‌های احتمال است. اگر 
| 
پیش‌بینی عملکرد پیش رو به طور قطع غیرقطعی اسست. اين بدان 
معنی ات که نتیجه هر گونه تحلیل مبتنی بر ورودی‌های نمایش 
داده شده با توزیع‌های احتمال» خود یک توزیع احتمال اسست. در 
محاسبه عدم قطعیت به روش شبیه‌سازی مونت کارلو بر خلاف 
روش‌های تحلیلی» اساس کار شبیه‌سازی پدیده مورد نظر به تعداد 
۰ بار می‌باشد. بدین منظور هر یک از پارامترهای ورودی 
(داده‌های تصادفی) به تعداد ۱۰۰۰ بار توسط روش‌های موجود تولید 
داده» شبیه‌سازی شده و با قرار دادن آن‌ها در مدل مورد نظر ۱۰۰۰ 
تکرار خروجی حاصل می‌گردد. سپس ضریب تغییرات پارامتر خروجی 
محاسبه و به‌عنوان معیار عدم قطعیت خروجی در نظر گرفته می‌شود. 
برای تولید هر پارامتر ابتدا بهترین توزیع آماری به هر یک از 
پارامترهای ورودی برازش داده می‌شود. به هر بار شبیه‌سازی» تحقق 
سیستم گفته می‌شود. برای هر تحقق, تمام پارامترهای غیرقطعی 
نمونه‌برداری خواهند شد که یعنی یک مقدار تصادفی از توزیع 
اختصاصی مربوط به هر پارامت, انتخاب می‌شود (۱۷). سپس, این 
سیستم با پارامترهای ورودی معین» در طول زمان, شبیه‌سازی خواهد 
شد. این شبیه‌سازی به گونه‌ای انجام می‌شود که در نهایت بتوانیم 
کارایی سیستم را محاسبه کنیم. این امر منتج به ایحاد تعداد زیادی 
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۰ نشریه آب و خاک. جلد ۳۵. شماره ۴. مهر - آبان ۱۴۰۰ 


نتیحه مستقل و جداگانه می‌شسود. که هر کدام بیانگر یک «آینده» 
احتمالی برای سیستم هستند. به عبارت دیگر یعنی یک مسیر احتمالی 
خواهيم داشت که سیستم با گذشت زمان آن را دنبال خواهد کرد. 
نتایج تحقق‌های مستقل سیستم به شکل توزیع‌های احتمالی 
خروجی‌های ممکن خواهد بود. در نتیجه» خروجی‌ها به صورت مقادیر 
تک نیستند. بلکه توزیع احتمال هستند. 


مدل‌سازی و نتایج و بحث 


با توجه به اینکه انتخاب مدل‌های ورودی در سامانه‌های هوشمند 
می‌تواند بر روی دقت جواب‌های حاصل از تحلیل تأثیر گذارده لذا سعی 
گردید پارامترهای مناسب و موّثر در مدل‌سازی بار رسوب معلق 
انتخاب گردد. جهت تخمین بار معلق رسوبی در ایستگاه اول از 
داده‌های روزانه مربوط به همان ایستگاه و با در نظر گرفتن ترکیبی از 
داده‌های دبی جریان و دبی بار معلق رسوبی با تأخیرهای زمانی 
مختلف و در ایستگاه‌های دوم از داده‌های مربوط به آن ایسستگاه و 


ایستگاه قبل استفاده شد. جدول ۱ مدل‌های تعریف‌شده در تحقیق را 
نشان می‌دهد. در این جدول 12۷ و 15 به ترتیب بیانگر دبی جریان و 
دبی بار معلق رسوبی می‌باشند و اندیس ) تأخیر زمانی را بر حسب روز 
نشان می‌دهد. همچنین در شکل ۲ نمودار مقادیر دبی جریان و رسوب 
معلق در دو ایستگاه انتخابی نشان داده شده است. لازم به ذکر است 
که ۷۵ درصد از ابتدای داده‌ها جهت آموزش مدل‌ها و ۲۵ درصد 
باقیمانده جهت آزمون مدل‌ها مورد استفاده قرار گرفت. 

جهت انتخاب تابع کرنل مناسب رگرسیون فرایند گاوسی و روش 
ماشین بردار پشتیبان مدل ۷15 با کرنل‌های مختلف مورد ارزیابی قرار 
گرفت. از معیار ارزیابی 115۳ جهت تعیین کارایی هر یک از 
کرنل‌ها استفاده گردید. نتایج به‌صورت شکل ۲ برای سری داده‌های 
آزمون ارائه گردیده است. مطابق با شکل ۲ مشاهده می‌گردد که در 
پی صا تترای کل تنم رل شساغ ماو یتفن ار 
معلق نسبت به توابع کرنل دیگر دارای خطای کمتری بوده و نتایج 
دقیق‌تری را اراکه کرده اسست. بنابراین» در تجزیه و تحلیل تمامی 
مدل‌ها از کرنل شعاع مداری استفاده شد. 


جدول ۱- مدل‌های تعریف شده در تحقیق 
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شکل ۲- نمودار سری زمانی دبی جریان و دبی رسوب معلق ایستگاه‌های منتخب 
عصمتاهای ۵60)60ز م6 10۲ مهد گنل )صمصصن0عی 060 یی 0صه 6عهه) ون 110 عمط ۵۶ عاماهر 1۳6 -2 عاع۲1 


سیدمهدی اقبیان. تخمین بار معلق رسوبی با استفاده از روش‌های هوشمند تلفیقی با در نظر گرفتن عدم قطعیت مدل ۰ ۴۸۱۲ 
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شکل ۳- پارامترهای آماری مدل ۳1۴ و ٩۷۷]‏ با توابع کرنل متفاوت 
امصم۱ امن ۱6 عام0مصهر 6۱۷۱۷۲ ره خ ما ۵۶ متمامصصه هن ام‌تاعتاهای 16 -3 مع۲1 
نتایج به‌دس تآمده برای مدل‌های میتنی بر کرنل خطا را داشته و مدل برتر است. مشاهده می‌شود که در حالت استفاده 
خر من ی زیمتو مایق ور اسان کی خریبان و بار عاق از پارامترهای دبی جریان و دبی رسوبی به‌عنوان تنها پارامتر ورودی 
رسوبی مربوط به روز قبل چند مدل تعریف گردید و کارایی روش‌های املل‌های ۷۱۱ و ۳۹ فایج دفعی ارانه شاه امشتت :ورهن‌های 
مبتنی بر کرنل مورد ارزیابی قرار گرفت و نتایج حاصله به‌صورت ترکیبی دو یا سه پارامتری دقیق‌تر می‌باشند. بین دو روش هوش‌مند 
جدول ۲ و شکل ۴ ارائه شد. با توجه به نتایج به‌دست آمده از سه معیار به‌کار رفته روش رگرسیون فرایند گاوسی نتایج دقیق‌تری را ارائه داده 
ازژیاین م | یرای سر خی آیسگاه مشب ناهد می‌ شود گذمدل است. اما در حالت کلی و با مقایسه نتایج مدل‌ها مشاهده می‌شود که 
5 با پارامترهای ورودی 95 ,۷ ,2۷ بیش‌ترین کارایی را روش‌های مبتنی بر کرنل از دقت مطلوبی جهت مدل‌س‌ازی دبی بار 
داراست. این مدل بیشترین ضریب همبستگی و تبیین و کمترین مقدار معلق به‌ویژه در نقاط ماکزیمم و مینیمم برخوردار نمی‌باشد. 


جدول ۲- نتایج ارزیابی مدل‌های مبتنی بر کرنل برای دو ایستگاه 
و ۷۵ 10۲ 0615مص 0عمدصا-مصعع۱ مامت اون 0۶ واابافع1 -2 12010 
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۲۳ نشریه آب و خاک. جلد ۳۵. شماره ۴. مهر - آبان ۱۴۰۰ 
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شکل ۴- نمودار مقادیر مشاهداتی و محاسباتی مدل 6۳18 برای دو ایستگاه (الف) مدلسازی با استفاده از داده‌های خود ایستگاه و (ب) بررسی 
ارتباط ایستگاه‌ها 


5ص 6 جدن 0ععوها مصتلم۵0صظ :(2) وعممتاهای 6۳۲۵ ۵۲ 06 مر خا ۳ ۵۲ مماه۲ ما60 20 ۷۵60 وهای عط) ۵۶ عامان 1 -4 عساهز۲ 
«تحاعصمتاعاه۱ مممتاهای عط) مصتامعنامه ۲ص خرط) مضه ماه 0 


نمودار مدل برتر برای سری داده‌های آموزش و آزمون برای هر 
دو ایستگاه با استفاده از روش 01۳ در شکل ۴ نشان داده شده است. 
همچنین جهت بررسی ارتباط بین ایستگاه‌هاء مقدار دبی رسوبی 
ایستگاه دوم بر اساس داده‌های ایستگاه اول به‌دست آمد و نتایج در 
انتهای جدول ۲ ارائه شد. در واقع مقدار بار رسوبی ایستگاه دوم برای 
حالتی که هیچ اطلاعاتی در مورد داده‌هایش در دسترس نیست بررسی 
گردید و از داده‌های ایستگاه اول تحت عنوان مدل‌های 2۷6 ۷۲7 و 
8 به‌عنوان ورودی استفاده سد. نتایج پیش‌بینی نان داد که 
درصورتی که از داده‌های دبی جریان و بار رس‌وبی به‌طور هم‌زمان 
به‌عنوان ورودی مدل استفاده شود نتیجه‌ی بهتری نسبت به حالتی که 
تنها از دبی رسوبی به‌عنوان ورودی استفاده شود به‌دست می‌آید. در 
این حالت مدل با پارامترهای ورودی 1251 ,152,12۷ عمکرد 
بهتری نسبت به دو مدل دیگر دارد. 


نتایج به‌دست آمده برای مدل تلفیقی -11۱۲0-6۳1۸- 0۱۷ 
اتقش 

در این قسمت. جهت افزايش دقت پیش‌بینی مدل‌های ٩۷۱۷۲‏ و 
15 از روش‌های پیش‌پردازش سری زمانی شامل تجزیه‌ی یکپارچه 
مد تجربی کامل و تبدیل موجک گسسته استفاده شد. بدین منظور» 


ابتدا سری‌های زمانی توسط روش تبدیل موجک تجزیه گردید. سپس 
با استفاده از روش ]۳۳۷ زیر سری جزئی حاصل از تجزیه 1۷۷۲ 
با فرکانس بالا دوباره تجزیه گردید. سپس انرژی تمامی زیر سری‌ها 
محاسبه شده و زیر سری‌های با انرژی بالا به‌عنوان ورودی در 
مدل‌های مبتنی بر کرنل استفاده شدند. اولین قدم برای تجزیه سیکنال 
مادر 104 با سه سطح تجزیه استفاده شد لذا دو سه سری جزئی و 
یک زیر سری تقریب حاصل شد. اساس روش ۲17۷۲۱ نیز تجزیه 
سیگنال به ۳/آهای مختلف و یک سیگنال باقیمانده است که در 
نهایت مجموعه‌ی این‌هاء همان سیگنال اصلی خواهد بود. روش تشکیل 
وقتی ادامه می‌یابد که تقریباً سیگنال باقیمانده ثابت شود. در شکل (۵- 
الف) زیر سری‌های حاصل از تجزیه بار رسوب معلق ایستگاه اول 
نشان داده شده است. مطابق با این شکل 1۷۳آهای دو و سه از 
زیرسری جزئی اول و 1۳های سه و چهار از زبرسری جزثی دوم» 
زیرسری‌های غالب حاصل از تجزیه با روش 2۷1 هستند. نتایج 
سری‌های زمانی با روش ۷۷۲( و تجزیه بیشتر آن‌ها با روش 


۱/0 


0037 
0۳040 
029 
)0(0(31 
0045 
)(048 
0032 
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۱۳ 


سیدمهدی اقبیان» تخمین بار معلق رسوبی با استفاده از روش‌های هوشمند تلفیقی با در نظر گرفتن عدم قطعیت مدل 


1 تا حدود زیادی باعث بهبود نتایج شده است. با توجه به نتایج 
به‌دست آمده از سه معیار ارزیابی مشاهده می‌شود که در این حالت نیز 
در ایستگاه اول مدل ۷5 با پارامترهای ورودی 1251 ,1۷۷۱ ,1(۷ 
بیشترین کارایی را دارا است. ملاحظه می‌شود که در حالت استفاده از 
پارامترهای دبی جریان و دبی رسوبی به‌عنوان تنها پارامتر ورودی نیز 
دقت نتایج تا حدودی افزايش يافته است. بااین‌حال» مدل‌های ترکیبی 
دو یا سه پارامتری دقیق‌تر می‌باشند. لازم به ذکر است که در استفاده 
از مدل‌های مختلف» عملکرد آن‌ها در تخمین مقادیر حداقل و حداکثر 
سری زمانی باید در نظر گرفته شود. مطابق با نتایج روش‌های مبتنی 
بر کرئل تلفیقی از دقت مطلوبی در مدل‌سازی نقاط ماکزیمم و مینیمم 
بار رسوبی معلق برخوردارند. همچنین در حالت تخمین دبی رسوبی 
ایستگاه دوم بر اساس داده‌های ایستگاه اول» نتایج پیش‌بینی نشان 
می‌دهد که مدل بر اساس هر دو داده‌های دبی جریان و دبی رسوبی 
مدل برتر بوده و با تجزیه داده‌ها کارایی مدل تا حدود زیادی افزایش 
یافته است. در حالت کلی» در صورت استفاده از روش‌های تلفیقی 
پیش‌پردازش دقت مدل‌ها تقریبا به میزان ۲۰ تا ۲۵ درصد نسبت به 
روش‌های منفرد هوشمند افزایش می‌پابد. 


۳۸۳ 


بررسی عدم قطعیت مدل 

برای بررسی عدم قطعیت. از دو فاکتور استفاده می‌شود که اولی 
درصد خروجی‌های مورد مطالعه است که در محدوده 957/۳۴۳ قرار 
داشته باشد و مورد بعدی میانگین فاصله بین باندهای بالا (260) و 
پایین (26) عدم قطعیت است (۱۲). در اين راستاء مدل در نظر گرفته 
شسه باند بارها الخرا شود [ هه پاررفر ای مظالهاد و احیال توزیع 
تجمعی مدل ها محاسبه شود. باندهای بالا و پایین به ترتیب ۰۲,۵ و 
۷,۵ احتصال توزیع تجمعی در نظر گرفته می‌شسوند. در سطح 
اطمینان مناسب باید دو نکته مهم در نظر گرفته شسود. اولا» 
باند 9590۳۳7 شامل بیشسترین داده‌های مشاهداتی باشد و دوماً 
میانگین فاصه بین باندهای بالا و پایین (1۳2610) از انحراف 
استاندارد داده‌های مشاهداتیکمتر باشد. می‌توان متوسط عرض باند 
فاصله اطمینان را به شرح زير محاسبه کرد که در آن 0 و 7 به 
ترتیب انحراف معیار داده‌ها و عرض متوسط باند اطمینان و ک تعداد 
داده‌های مشاهداتی است. 


جدول ۳- نتایج ارزیابی مدل‌های مبتنی بر کرنل با پیش پردازش داده‌ها برای دو ایستگاه 
کصمناهاه ۳۲۵ ۵۲ مصمی ۵۵اه )02 ۱۲۱6۵ واع0 مج 0ععوها-امصم عصتامتاوی ۵ وااناعع1 -3 1016 
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(82) 116 زمان (روز) 
شکل ۵- (الف) زیرسری‌های حاصل از تجزیه بار رسوب معلق ایستگاه اول. (ب) نمودار مقادیر مشاهداتی و محاسباتی مدل تلفیقی 2۳1) برای دو 
ایستگاه با استفاده از داده‌های خود ایستگاه و (ج) بررسی ار تباط ایستگاه‌ها 


۵ 20 6۳۷60وآن 6 ۵۲ واماص :(0) با ممتاهای عط) ما ممزم‌وهامای 0ععممجصمعع0 0هم۱ )صمصت0عی 0260هص۵صرونای 16 <(ه) -5 ۲1۵0۲۵ 
متصاعصمت)هاه۲ عصمتاهای عمط عصتاهعتامه 10۲ 0(۰) 200 مصاجل هه و صمتاهای عط) عصمه فصمتاهای ۳۲۵ 0۲ 0۵ص ۳ 0عا عم ۵ دم(۲۵ 


سیدمهدی اقبیان. تخمین بار معلق رسوبی با استفاده از روش‌های هوشمند تلفیقی با در نظر گرفتن عدم قطعیت مدل ۰ ۴۸۵ 


جدول ۴- نتایج تحلیل عدم اطمینان مدل‌های 618۳ و ۲۷۷۲-۳۷۲-6۳1۴ 
)۱6 ۳۸)-(۲/۲,/۷۲1- ۲۷۷ 20 6۳۲۸ 66 ۸0۲ عناوم وزووافه زاصنه نمض -4 12016 


معیارهای ارزیابی دهانت ۳۵۲00۵۵۵ روش 
0-0۳ 9۳۳1 / ۰ ۰ 95۳۳۲ ۲ ۰ ۰ 95۳۳۲ 00( 
025 1-140 ۵ 0242 01-0 0341 12.10 0۳1 
3 3 
گ ک ک 
0۰110 ی 5 0105 ۵( 0093 93250 موم رصن 
ام 
7۳ ش جدید جهت بهبود دقت پیش‌بینی با ق پرداخته شد. 
)6 ۱ روش جدید جهت بهبو پیش‌بینی بار رسوب معلق پر 
0 مدل‌های ورودی متفاوتی با استفاده از داده‌های مربوط به دبی جریان 


۱۰ > > ,)6 99۳۳۱ ره نعه 3:۵ و دبی رسوب و با در نظر گرفتن دو نوع مدل‌سازی با استفاده از 
داده‌های مربوط به خود ایستگاه و داده‌های مربوط به ایستگاه قبل 


نتایج بدست آمده برای تجزیه و تحلیل عدم قطعیت در جدول ۴ 
تعریف گردید و مورد بررسی قرار گرفت. نتایج حاصله نشان داد که 


نشان داده شده است. براساس مقادیر بدست آمده برای شاخص‌های 
۲ماع0۳ و [9590۳۳1 مشاهده می‌شود که بیشتر مقادیر تجزیه سری‌های زمانی با هر روش ۷۷۲ و تجزیه بیشتر زیرسری‌های 
ی 
همچنین» مقدار شاخص 1-۳26107 برای مجموعه داده‌های آموزش و می‌شود. مطابق با نتایج حاصله. برای مدل‌های تعریف‌شده براساس 


آزمون کمتر از انحراف استاندارد داده‌های مشاهداتی است. بنابراین» 
براساس نتایج می‌توان بیان نمود که مدل‌سازی بار رسوب معلق با : 
استفاده از مدل تلفیقی 06۳۴ دارای درجه عدم قطعیت محاز می‌باشد. ورودی 125۱ ,۷۷۱ آری,5] با داشستن بپیشتر > و 10 و کمترین 

مقادیر خطا مدل برتر به‌دست آمدند. مشاهده شد که در صورت 


استفاده از هر دو روش‌های تلفیقی نتایج دقیق‌تری در تخمین بار 


داده‌های خود ایستگاه» مدل ۷۲5 با پارامترهای ورودی ,1۱۷7 ,1۷7 
51 و در حالت بررسی ارتباط بین ایستگاه‌ها مدل ۷18 با پارامترهای 


۱ رسوب معلق حاصل شده و دقت مدل‌ها حدود ۰ تا ۲۵ افزایش 

پیش‌بینی بار معلق رسوبی به دلیل اثرات آن روی شاخص‌های می‌یابد. همچنین نتایج عدم قطعیت نشان داد که روش تلفیقی -۲۷ 

کیفی آب تقلیل گنجایش مخازن و تغییر در مورفولوژی رودخانه‌ها از 1۲:۷۴ دارای درجه اطمینان مناسبی در مدل‌سازی بار معلق 
اهمیت ویژه‌ای برخوردار است. لذاه در این تحقیق به پیش‌بینی بار رسوبی است. 


رسوبی معلق با استفاده از روش‌های تلفیقی پیش‌پردازش به‌عنوان 
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,۷6۵۲5 ۲۵۵۵۵۲ 1۳ رصق تعطان مطا ۵ .همه تمطونط 16 05مطاعصه تما موبا م6 26665527 ور )1 تناو متأنیجل بط 
۵ 001۲۵۱6۴ 1۳0108216 20 متاباهتل بط مط مصتاهعتاوع7 1۳ ما 0عناممره موه 2۷۵ فمطمه2۵020 9001 ه)۱6 عطا 
ما عصا مرحم 61160017 و9 ما ماو صفعها وقاچ 7۵« ممتلومممصممعل اقصونه عمز 1270۱ واممصه ۲۱۱۷/۵۱۵۵ 
6۵ ۱1۳۱624 ماطاصهعص نجتممصمام‌همن ما 04162060ص قق. رولمطاعصه عمتاعتهي فعتیعو مصصتا 0۶ «مهتنا0ع2 
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.1( 
امه ۵1 11601 مطا 0ص ممصمصمصمطم ممتومیته مه اممحصتلعد مطا که همم فطا مه فتل رتلتتاد فتظطا 1 
مج (6۱۷) مصتطممحه ۷۵۵۵۲ آ0۲مجنای طا عصمله وعمط)فمصظ فصلوومع۲۵م-0۲۵ وملنعد محصلا . رعصتامصتاوه صا فتمامصطه221 
۵ ۱020 ]ممصصنلع6ز 060ممم‌ویاو متقحتادع 60 ۱۵۹۵0 ۱۵۲۵ وعظ)ع20۵۲02 0عقهه امصریع1 رعرطن) ممتووع۲۵۵۲ ووعع0۲۵0 صهلفوا02) 
۶ 02960 0011060 ۷۷۵۲۵ دام )صعععن رعوممتنام فنطا ۲۵۲ .عفصمتاهاه متامصم وه عتانا۵عفممم ۵ 26 ۲۱۷۵۲ تاه 
005 20 عصلووع0۲۵6-۲06 0ماهوماصا گه دوه ما ره0۲ع۷۵۲ .متا وتمامه فده ومامتانهج اممصصتلعو م24 متلنیهتل بط 
۶ 36 ات۳۲ 4ماهعلاوه12۷ ۵۵ مممتاماومهع0۵0۷ مه حملاهافماصا ۵1 وعاقاه ۸۵ طا ممطامقظ. هصلووعع2۳۵ 
٩0160160 40‏ ۷۷/۵۲۵ و6۵۲1۵و-او ۳60۱2۵6 طونط مطا رصقطا عصتووععم1م-016 0202 ۲۵۲ ۵0ویا ۷۷25 00ظ)عح (۱۷۷۲) نومه[ 
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مایم ۱۷۲۱۵۵۱۵ ,م۵0 عمتنمم‌داه عط ۶ باتازطدتام۲ مط) ووعوفه مه رصمتت20 ما .فاعم 0عقع‌-1عهزم۱ ۶۵۶ وانامرصاً فج 1۳200560 
ص راصتماهعصه ۵1 062۳66 ماطمتنوع 2 20 ۱۵061 کظ0۳) مطا )ها 0760و فاناوع م1 .۵0و 2۵ فلوراممه اصتماهعص 
018۰( 
ععصهادنل م1 .۵0و وه ۳۷۵۲ متمماهوتام۲۱ ۵۲ فحمتافاه ۳۵ ]۵ عاجل یناه عتط ح1. :فل‌مطاع]۷ 0ص2 واهزه) ۱۷2 
فص ,مااهکنطهووه]۱۷ رصماطعتدظ کمعین عدعه 1۱062160 و1 ممتاه)ه افتت مط1 بح 50 «امتمصتدمتم2 قه فصمتاهاه طمع۵ 
6 ۲650۵6۵0۷۵۱۰ روعانص معمتانو 634 20 282 15 معمتاهاه عطا 10۶ معته صلجوه مظ1 ,عامتامعصجمن طا وا ممتاهاو ۹660840 عط 
کقطا ممتام‌حصناووح مط مق 0عقهه مه امعم کزظن) 0ص 5۱۷۱۷۲ .صمتاهای 966080 فطا م6 صمتاقاه اف فطل ص۳0 و1 هم 1100 
٩0601818 2‏ 0۴ ۷۷۵۷ 2 276 0۲06689565 0209912 .تعطان طمهه )باهطاه صمتهصمگصا موم 10بامطای فصمتاهتعوهاه عصعمهر20 
0ص صفمصه مدمه ممتاطاتتافت صوتوفنیجی فط گم ممل)ه2نهتعصمع لهعنهه م و1 فنط1. .عمهمو ممتاعصظ هن امعت تمتم 
امطمتمصط 20 0242 فط عتامطاج معلع(مص۱ منم ما مان .۲۵8۵۵00 بکتتاقمه 24 ۷۵۵۵۲ 2 276 00۷2112800 
ما وعلصه م) 01اه یه و01مع ممتووع1عع۲ ظن) 4صه رحمتا2112 تمصع ۶0۲ ۲60۵1۲۵0 و1 0۴06695 ۷۵11020109 مظ رقعآممع0 ه۵۵ع0 
۷ ۱۵:19 2 وقعوا (۷۲۲۲) صرمقصه 1 ]ما۷۷۵۷۵ اتامرصا اوما معط مه عصتل‌جممعمتمی صمتااماتتاوتل و۷تامت0ع۳۳ مطا 
عصل20001 مصصنا کع۷۲ن 00عصعدان مها هم صمتامصن مه مامان 11۵ مط 1 ,ممتوومممم اقصوه صا (۱۷2۷۵1۵6 تمطامع) صمنامعصن] 
-عع127) فممتمصتد۵۲0م2 020 مصا موممررمعع 1۱ تمصرعنه مط) ,۷۷۲ عصلمیا معا .وومصامهم‌مههن 200 عمقمه اقصعنه مطا ما 
0 01۳000960 ۷۷/25 ۱۳۱۱۷۲۲۱ .واممممموصمم قمم‌ممم‌صومع ملهعو-211ظ2) 06021160 4ص قصمممم‌مم مممنمم- 1۵0۳ ۵۲ ملهعو 
صومصط مطا وج ۲۳۲۳ مبایا عمط وعتمعمه طمنط (۱۲ظ) ممتاتومم‌معع0 معط لهمزت1م‌صصن ۵ مناقعا لدنص مق فص 801۷۵ 
فنص ۶ وعامط مانط/ 2 فبتام اقصعته مطا ۵۶ عناوم ممت)تومم‌جطمعع0 فص گم فاعتعومی مت مق علمت ۵۶ قاماجمعصه صه ۵۶۴ 
,۲691006 2 20 مومتاعصب مل‌مصط متفصتاصا ماتص ملصا تقصعه عآمصمن ره 0۵601۴9۵098 مه ۵0فنا ۵8 هی ۱۸۲۲۸ .ع4تالا2۳0 
40 2021۷776 ۲0 موه مه اقا معمتهاه‌م ممممتاو 20 فامه‌ممم‌مومع «مممتامع۲ موه طز فافم‌طاناه طا عصتالناوع۲ 
۱01886 10۲ ۱۷560 06 620 11 2 و ممتالومم‌جطمعع0 ۵ ام0مصه مزونه ۵۶ معنقجع؟ اصما۲مصا تمطامصن بافوعع۵۲] 
نع معط مز .عممتامن0عزن (م0مصظ ۵ «مجتتامعه مطا عصامتمصا مد متامعه وه هم من روعترعو مصصتا م0 ۵0 
عم ح م1 .اصلقانهعصنا واملموصه مص 22172 م6 صرح وومصامیاطام عمط )وم م6 12860 27 واصمصمام ۳۵ رل‌مطامصظط ولو 2821 
6۵ 2۷6۵۲۵826 6ظ) 15 عمه ۱۵2 مطا همه 95۳۴۳۴۱۲ ۵۶ موصده: عمط صد میم )نها فتاماتان 0160و عطا ۵۶ ععفادهعتعم معط و1 
۷ 1۲۱۱ 08 یباهو [06م 0مع10فومع فطع رتقعع منط مر .عصهه )معا ریک و10 هه ریک تممرمرنا مط) معع۵ 


صو رتدطه نمهب دنا ۸2240 متصصهاو1 رطمصحظ فطم رعمتهعمتعمظ از ۵۶ امعصتدمعن1 -1 
(۵121-2027.26.11) مهاطامطمه فد :اتف میمصت طمصومتزمن) -۶) 
1 3 :]1۵ 


۸ نشریه آب و خاک. جلد ۳۵. شماره ۴. مهر - آبان ۱۴۰۰ 


۰ 08 0615مهص مط) گم بوناتهدهامت ممتتاطاتتاوتل 1260۷۵ نیم لهعتنژم‌جصم فطع مه ,رنه فنطا صر وعحصت 1000) وفحصت 
بطمتاطاتافتل وبتاقاصته فطع ۵۶ ممتاتازطه‌طامتم 97.52۵ 280 2.560 106۲60قجهع عتع فده م1 صه. تممرهرنا. فظ 1 
0۰( 

مصتصماعل 0صد2 ۳۵0۵15 0عاوع) عطا 0۶ ممصمممتتهج مصا مزب مج متقتالهبه ما م۵0 1 خومتکفیه کزظ 0صد و6اناعع1 
ممتامصتصصهعاه 1 بر اصمه‌للهم ممتاهاهترمی مه همه ممصمطمگهم معط رولع0مصظ 96160160 عط ۵۶ «مجساعع2 عطا 
مطا )2 10162060 فااناوع 0مصتق‌طاه مط1 ,لمونا موب (رظا(فع) وم ممتاوو مصقع] ۴۵0۲ هه ,۲0 )ممملهم) 
۸ 15۵ ۲۳6 .۵06015 ۳۲۴ مج ۱۷۲ 5۱۷ 16مصده عط) مقطا ععطمنط وه واع0مهه 72060عمص1 0عناض(2 عطا ۵۶ رمهتنا0ع2 
وم لو همه اصتدانهعص مطا 10۶ دالناوع؟ 0معصتق)0ا۵ 16 ۹۵۰ 25 6۵ 20 همع07عه ماه مه فطا 060762860 و8مطامعظ 
وم ها عصح‌ها ۳۳۲ 95 عمط منز ۷۵۲۵ ومبله۷ ۵۲6010)60 20 ۷60تعوهاه فطع عصناعمصظ 1020 )ممحصعو 060جهم‌فیای صا فا 
عط صقطا تعالمصه ۷۷۵۲۵ ماعع02128 )ومع 0ص2 صتق ۶۵۲ 1-1۳2000۲5 ۵۶ اصتامصصيح مه ق 0صنام م۷2 11 رته ۷۲۵۲60۵۷ .وعووع مطا 0۶ 
0 1386 )12 10660 9 نامع 16 رقتلتافعع معط من 0و2 ,10۵۲۵]01۵ ,0242 ۲۷۵0عوما۵ فطل ۵۶ 0071200 0270صماه 
۰ ۶ 162166 2110772016 مج 60 160 8۵061 > 0۳-(1-۳:۳/۷۲۲ ۷۷ 1006212060 ۷1۵ م0061 اصمحصتل86 

٩۱۷ ۱۷ 5‏ 4ج 0۳ 212000ع0ض1 اقطا ۲۵۷۵2۱60 بمجتتامعج "داعم 00۷7010۵۵0 مطا که مموتتهممصمن خصمتفسمومن) 
صمصصنلعو 060همم‌ویاه فط) نامهتم ها وام0مصه 6۷۱۷۲ مه عاطن ملمصته اه همهم عمصممصزمتوم تعطونط 24ظ 
018 .0 25 0 20 ۵6076۵ 6۲1۵ میم معط 060۲62800 «امام0۲0 ۵۵ 05مطصامقطر و عوفطا ۵۶ عویا فظ 1 ,عععقطمول 
۵۹ / 1۳0۷ عط زب 9001 مصا بلق متفه ما چم 0عفهه 06۷6010۵۵0 م۱۷۵۲ )2 وام۵0هظ مطا 10۲ رقلنافع؟ فطا ما 
عط ۳/16 9۵06 معط رقصمتاقاه مط ممعباه متامممتاهام عمط عصتاهعلاوع1۵۷ ۵۶ مفهه مطا ما 240 تون 20 رن نا 0 
معط رولوررلممه «اصتهعهعص مط جم 0عفهه ,وله .عاممهص عملوهناد م۷۵۲۵ رون 200 ونع 1۳ رم‌ونا 0 اهتدم )1۳0 
0 1 ,۳۱۵۱۷۵۷۵۲ .عصت1ع0 معط )مرمحصصیلعو 60ممموتای ما راصتدانهعصن ۵1 ع7عع0 2110772016 مه 20ظ ۵061 ک0۳) 1۳16812060 
فتطا ۵۶ تاه وععص12 0242 عصلونه دعتلنتاه ماع ر۵ع۵۲۵10ظ۲۳. .وام0مه ۷۵ الفصمو 0242 مه 05مطاممطر ۱960 مط عقطا 0مامظ هط 
اصمصصلمو 060مممرمیاه متقصصتاوع م6 ولععمصظ معط گم مترمصه ما مصتصصعاع مه نام 2160 ۵۵ 10تامطو هه 1610 مه 0و 
۰ 0 وصمنایجمع ۲62 عط) مد 1020 


0 ,۲6-۵۲00۵6551172ظ ر1681695101 2۳068۵85 299120 رطملالو0م‌جرمعع. 8مظ له)صمصصهصد 1۵۲۱۷۵۲8۰ 
اصعصصتل6و 


